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Jummary 

A  comparison  study  was  performed  on  four  gear  dynamic 
nalysis  computer  programs  developed  under  NASA/Army 
ponsorship.  These  programs  are  GRDYNMULT  (a 
nultimesh  program  applicable  to  a  number  of  epicvelic 
\  stems),  TELSGE  (a  single  mesh  program).  PGT  (a 
nultimesh  program  applicable  to  a  planetary  system  with  three 
ilanets).  and  DANST  <a  sing!*1  ;,.^>h  program).  The 
apabilities  and  features,  input  and  output  options,  and 
echnieal  aspects  of  the  programs  were  reviewed  and 
ompared.  Results  are  presented  in  a  concise  tabular  form, 
’arametric  studies  of  the  program  models  were  performed  to 
nvestigate  the  predicted  results  of  the  programs  as  input 
tarameters  such  as  speed,  torque,  and  mesh  damping  were 
■  aried. 

In  general,  the  program  models  predicted  similar  dynamic 
oad  and  stress  levels  as  operating  conditions  were  varied.'- 
-lash  temperature  predictions  from  programs  GRDYNMULT 
ind  TELSGE  indicated  similar  trends:  howe'  er,  actual  values 
■vere  not  in  close  agreement.  The  program  GRDYNMULT 
vas  found  to  be  the  most  versatile  in  system  size.  type,  and 
tnalysis  capabilities.  The  programs  DANST.  TELSGE.  and 
3GT  are  more  specialized  for  specific  systems;  however,  in 
ipeciftc  areas  they  provide  a  more  detailed  treatment  than 
3RDYNMULT. 

Introduction 

Since  the  late  19th  century,  gearing  has  become  the  simplest 
ind  most  efficient  means  of  transmitting  mechanical  power. 
Hears  can  be  found  in  almost  every  application  involving 
nechanical  power  transfer,  and  arc  usually  considered  a 
.ritical  link  in  the  power  chain  of  that  system.  Because  of  this, 
tear  designers  are  highly  concerned  with  gear  I  if*’  -ind  rcli- 
ibility.  In  industrial  applications  this  concern  is  alleviated  to 
iome  degree  by  over  designing  the  gears,  sacrificing  cost,  and 
ncreasing  weight.  However,  in  aerospace  applications,  where 
veight  and  size  are  premiums,  gear  systems  are  usually 
lesigned  close  to  their  projected  limits.  As  a  result,  a  number 
if  computer  programs  have  been  developed  in  an  effort  to 
iredict  parameters  such  as  dynamic  load,  surface  damage,  and 
.urface  temperature,  that  are  integral  factors  in  various  gear 
ailure  modes.  Several  of  these  pro--  mis  have  been  developed 


through  NASA  Lewis  Research  Center  under  NASA/Army 
sponsorship. 

Of  all  the  gear  dynamic  programs  developed  at  NASA,  the 
programs  TELSGE.  GRDYNMULT.  PGT.  and  DANST  are 
the  most  widely  used.  TELSGE  was  developed  to  study  the 
effects  of  input  parameters  such  as  speed,  load,  and  lubricant 
oil  type  on  predicted  quantities  such  as  dynamic  tooth  mesh 
loads,  surface  temperatures,  and  lubricant  film  thickness  in 
a  single  mesh  system  (refs.  1  and  2).  Gear  failure  modes  such 
as  scoring,  pitting,  and  lubrication  failures  are  directly  related 
to  these  predicted  parameters.  GRDYNMULT  was  developed 
to  predict  parameters  such  as  tooth  mesh  loads,  tooth  stresses, 
and  surface  damage  factors  under  a  variety  of  input  conditions 
for  a  single  mesh,  or  multiple  mesh  cpicyclic  system  (refs.  3 
to  5).  These  parameters  have  a  direct  effect  on  failure  modes 
such  as  tooth  breakage,  scoring,  and  pitting.  The  program  PGT 
was  developed  for  the  dynamic  analysis  of  a  three  planet 
planetary  gear  system  under  a  variety  of  input  conditions 
(ref.  6).  The  magnitude  of  the  dynamic  mesh  load  output  from 
PGT  indirectly  influences  the  probability  of  tooth  failure  by 
breakage.  The  program  DANST  was  developed  to  study  the 
effects  of  input  parameters  such  as  tooth  profile  modifications 
and  external  shaft  and  mass  magnitudes  on  predicted  dynamic 
loads  and  stresses  ot  a  singie  mesh  system  (rets.  7  to  9).  The 
tooth  root  stress  parameter  predicted  is  a  critical  factor  in 
determining  gear  failure  through  tooth  breakage. 

The  purpose  of  this  study  is  to  provide  a  comprehensive  guide 
on  the  capabilities  and  nature  of  results  obtainable  from  the 
four  gear  dynamic  programs  introduced  above,  and  to  provide 
some  program  verifications  through  direct  comparisons.  The 
report  is  divided  into  two  main  sections.  The  first  section 
reviews  the  capabilities,  input  and  output  options,  and  technical 
aspects  of  the  programs  studied,  and  presents  the  results  in 
a  concise  tabular  form.  The  second  section  reviews  comparison 
runs  that  were  performed  to  compare  the  results  obtained  from 
°ach  program  using  common  input  models  and  parameters. 
Finally,  some  concluding  remarks  are  presented  which 
generalize  the  results  of  the  total  comparison  study. 

Program  Features  and  Models 

Research  on  each  program  was  conducted  to  obtain  the 
general  and  technical  fetures  of  the  programs  on  an  individual 
and  collective  comparison  basis.  Program  features. 
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capabilities,  and  options  were  tabulated  in  an  effort  to  provide 
an  easily  accessible  reference  base  for  potential  program  users 
Table  I  presents  some  general  information  on  each  program 
such  as  system  sizes  and  types,  gear  types,  and  supporting 
documentation.  Table  II  gives  a  direct  comparison  among  the 
programs  of  the  type  and  nature  of  the  parameters  calculated 
by  each.  A  comparison  of  the  input  options  available  for  each 
program  with  some  basic  descriptions  of  these  options  are 
presented  in  table  III.  Finally,  table  IV  gives  information  on 
the  printed  and  plotted  output  options  available  with  these 
programs.  In  the  following  sections  general  program  features, 
as  presented  in  tables  I  to  IV.  are  discussed,  along  with  the 
various  analytical  models  used  in  the  programs. 

General  Capabilities,  and  Features 

Program  PGT.— The  program  PGT  (dynamics  of  Planetars 
Gear  Trains)  (ref.  6)  is  a  gear  dynamic  analysis  program  for 
a  three  planet  planetary  spur  gear  system.  PGT  is  capable  of 
modeling  a  planetary  gear  train  with  input  and  output  shafts 
and  masses.  It  calculates  dynamic  mesh  loads  and  combined 
stiffness  for  each  mesh  as  a  function  of  roll  angle.  PGT  also 
calculates  the  sun  center  movement  in  the  plane  perpendicular 
to  the  sun  gear  axis.  Along  with  the  standard  input  parameters, 
such  as  tooth  geometry,  torque,  and  speed,  other  parameters 
can  be  input,  such  as  profile  errors,  sun  center  stiffness  and 
damping,  etc.,  as  indicated  in  table  III.  The  major  features 
of  this  program  are  its  ability  to  include  input  and  output 
peripherals  in  the  analysis  and  to  calculate  the  movement  of 
the  sun  gear  center.  The  major  limitation  of  this  program  is 
that  it  can  only  be  applied  to  a  three  planet  system.  Sample 
plotted  outputs  of  PGT  are  given  in  figure  1 .  The  first  two 
plots  represent  the  dynamic  load  factor  for  the  sun/planet  and 
ring/planet  mesh  associated  w  ith  planet  number  1  of  run  1  in 
table  V.  The  dynamic  li  ad  factoi  represents  the  ratio  of 
dynamic  to  static  tooth  load,  and  is  commonly  used  when 
plotting  dynamic  mesh  loads.  The  sun  center  movement  plot 
is  the  actual  displacement  of  the  sun  center  through  one 
complete  steady  state  revolution.  It  should  be  noted  at  this  time 
that  program  PGT  is  not  in  an  easily  runable  format.  Some 
work  would  be  required  to  revise  the  program  to  a  more 
standard,  commercially  acceptable  status. 

Program  GRDYNMULT.  —The  program  GRDYNMULT 
(Epicyclic  Gear  Dynamic  Analysis  Program)  (ref.  3)  is  a 
dvnamic  analysis  program  with  the  capabilities  to  model  a 
variety  of  gear  types  and  gear  train  systems.  GRDYNMULT 
is  capable  of  modeling  single  mesh,  planetary,  star,  and 
differential  systems  with  a  maximum  of  20  planets.  This 
program  can  model  spur  or  helical  gear  types,  along  with 
involute  or  buttress  tooth  forms.  GRDYNMULT  is  capable 
of  calculating  a  number  of  variables  such  as  dvnamic  mesh 
load>.  tooth  root  stresses,  hertz  stresses,  flash  tempeiatures. 
etc.,  as  shown  in  table  11.  As  illustrated  in  table  III. 
nonstandard  parameters  such  as  tooth  spacing  errors,  tooth 
profile  modifications,  sun  center  stiffness  and  damping,  etc.. 
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fb)  Ring'planct  tooth  loads  mesh  1. 
fc)  Sun  center  movement. 


Figure  I  — PGT  sample  plotted  output  Input  torque  =  33.0  N*m;  input 
speed  -  4000  rprn. 


v  input  in  the  program.  The  ma  jor  feature  of  this  program 
»  variety  in  the  type  of  calculations  available,  and  the 
ber  of  gear  train  systems  it  can  be  applied  to.  The  major 
:ation  of  this  program  is  that  it  cannot  include,  in  the 
imic  analysis,  me  effects  of  input  and  output  peripherals 
.'ally  present  in  actual  gear  systems, 
imple  output  plots  from  GRDYNMULT  are  given  in 
res  2  and  3.  These  plots  are  for  the  ring/planet,  sun/planet 
h  associated  with  planet  number  1  of  comparison  run  1 
ible  V.  The  first  plot  in  each  figure  is  the  dynamic  load 
ir  for  the  mesh.  The  PV  plot  represents  the  product  of 
ocal  contact  stress  and  the  sliding  velocity.  The  PV  product 
sed  in  analyzing  surface  damage  possibilities,  such  as 


scoring.  The  Hash  temperature  plots  represent  the 
instantaneous  gear  surface  temperature,  and  the  hertz  stress 
is  the  local  contact  pressure.  The  planet,  ring,  and  sun  gear 
stress  plots  refer  to  the  tooth  root  stresses.  The  plots  associated 
with  GRDYNMULT  appear  different  from  those  of  other 
programs  because  GRDYNMULT  presents  only  half  of  the 
tooth  contact  cycle,  and  the  plot  includes  more  than  one  tooth 
pair  if  more  than  one  pair  are  in  contact.  Subsequent  plots  from 
GRDYNMULT  have  been  replotted  for  easier  comparison 
with  the  other  programs. 

Program  TELSGE. —  The  program  TELSGE  (Thermal 
Elasto-Hydrodynamic  Lubrication  of  Spur  Gears)  (ref.  I )  is 
a  dynamic  analysis  program  for  a  single  mesh  spur  gear 
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Figure  2  —  CiRDYNMl’LT  sample  plotted  output  of  sun  planet  mesh  1.  Input  torque  r  IV1)  N«m.  input  speed  -  4(MM)  rpm. 
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fa  I  Dvnamic  load  plot.  fh)  PV  plot. 

U)  Contact  stress  plot  ( d )  Flash  temperature  plot 

(e i  Ring  gear  tooth  mot  stress  plot.  ft)  Planet  gear  tooth  root  stress  plot. 


Figure  3.  —  ORDYNMl'I.T  sample  plotted  output  ol  ring  planet  mesh  1  Input  torque  =  33. N*m;  input  speed  =  4(HK)  rpm. 


system.  As  illustrated  in  table  II.  TELSGE  is  capable  of 
calculating  variables  such  as  film  thickness,  flash  and 
equilibrium  surface  temperatures,  dynamic  mesh  loads,  and 
hertz  stresses,  etc.,  which  are  important  parameters  in  gear 
tooth  surface  failure  models.  TELSGE  predicts  fatigue  life  of 
the  gears  based  on  these  calculated  variables.  Additional  input 
parameters  for  IELSCF  include  tooth  profile 
error  modification  array,  thermal  and  viscous  properties  of 
the  lubricant,  etc.,  as  seen  in  table  III.  The  major  feature  of 
this  program  is  its  comprehensive  treatment  of  the  dynamic 
and  thermal  effects  of  the  lubricant  on  the  resulting  life  of  the 


gears.  The  major  limitation  of  this  program  is  that  it  applies 
only  to  a  single  mesh  system.  Sample  plotted  outputs  of 
TELSGE  are  given  in  figures  4  and  5.  The  plots  shown  were 
constructed  using  a  postprocessing  graphics  program,  as  the 
current  CRAY  version  of  TELSGE  does  not  have  a  plotting 
routine. 

Program  DANST.— The  program  DANST  (Dynamic 
ANalysis  of  Spur  gear  Transmissions)  (ref.  9)  is  a  dynamic 
analysis  program  for  a  single  mesh  spur  gear  system.  DANST 
is  capable  of  modeling  a  system  with  input  and  output 
peripherals  included  in  the  analysis.  As  illustrated  in  table  II. 
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(a)  Dsnamic  load  plot, 
it)  Combined  tooth  stiffness  plot. 

Figure  4  —  TF.LSGH  sample  plotted  output.  Input  torque  =  203.4  N«m: 


tb)  Contact  stress  plot, 
t d )  Film  thickness  plot. 

input  speed  =  6000  rpm. 


(a)  Pinion  surface  temperature  plot.  (hi  Gear  surface  temperature  plot 

let  Pinion  flash  temperature  plot.  id)  Gear  (lash  temperature  plot. 


Figure  5  — TFLSGF  sample  plotted  output  Input  torque  -  203.4  N*m:  input  speed  =  6000  rpm. 


\NST  is  capable  of  calculating  dynamic  mesh  loads,  root 
esses,  combined  stiffness,  etc.,  as  a  function  of  contact 
sition.  Along  with  standard  input  parameters.  DANST 
ows  input  of  a  user  defined  tooth  profile  deviation  array, 
tndardi/ed  tixith  profile  modifications,  input  and  output  shaft 
d  mass  data,  etc.,  as  seen  in  table  III.  The  major  feature 
this  program  is  the  detailed  tooth  profile  error  modification 
put  available  to  the  user.  A  major  limitation  of  this  program 
that  it  applies  only  to  a  single  mesh  system.  Sample  plotted 
upuis  from  DANST  are  given  in  figures  6  and  7.  As  seen 
ligure  6.  DANST  provides  a  plot  of  the  Fourier  transform 


of  both  the  static  transmission  error  and  the  dynamic  tooth 
loads.  These  plots  can  be  useful  when  comparing  the  analytical 
results  with  test  results  in  the  frequency  domain. 

Program  Models 

Dynamic  models.— Vo  describe  the  dy  namics  of  the  systems, 
each  program  uses  differential  equations  of  motion  based  on 
mathematical  models  simulating  the  various  masses,  springs, 
and  damping  present  in  the  actual  systems.  The  mathematical 
model  used  in  PGT  is  shown  in  figure  X.  As  depicted  in  this 
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(a)  Stalk'  and  dynamic  toad  plot.  <hi  Tooth  root  stress  plot. 

Ki  Humor  transform  ot  static  transmission  error  id)  Kuiner  translorm  «>l  dynamic  tooth  loads 

Hiinre  ft  DANST  sample  plotted  output  Input  torque  20.V4  N*m.  input  speed  2000  rpm. 
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figure.  each  mesh  is  represented  by  an  equivalent  spring  anil 
dashpot.  The  spring  represents  the  combined  silliness  of  the 
gear  teeth  in  mesh,  anil  the  dashpot  represents  the  resulting 
mesh  damping  The  springs  and  dashpots  shown  at  the  sun 
eenter  are  present  to  model  the  flexibility  and  damping  oi  the 
sun  gear  shaft  and  hearings.  The  stiffness,  masses,  and 
damping  associated  with  the  input  shaft  and  driver  and  output 
shaft  and  driven  device  are  also  included  in  the  model.  Figure 

9  illustrates  the  model  used  in  the  program  GRDYNMF'FT. 
The  mesh  stiffness  and  damping,  and  sun  center  stiffness  and 
damping,  are  presented  similarly  as  in  the  PGT  model  As 
seen  in  figure  R.  additional  springs  representing  flexibilities 
between  ring  gear  rim  segments  and  between  planet  carrier 
segments  are  included  in  the  GRDYNMIT.T  model.  Figures 

10  and  I  1  represent  the  models  lor  programs  TF.l.SGF  and 
DASG.  respectively.  As  seen  in  figure  I  I .  the  DANST  model 
includes  the  mass  and  elastic  data  of  the  input  and  output 
peripherals.  Again,  the  mesh  springs  represent  the  combined 
stiffness  ot  the  gear  teeth  in  mesh.  For  a  more  thorough 
description  of  the  individual  models  and  the  iterafvc  methods 
used  to  solve  the  resulting  differential  equations,  refer  to  the 
supporting  documentation  for  each  program  as  given  in  table  I 

Tooth  stiffness  models.  —  To  model  the  complex  stiffness 
ot  gear  teeth  during  mesh,  all  of  the  programs  use  a  nonlinear 
tooth  compliance  model.  Programs  TFFSGF.  GRDYNMl  FT. 
and  DANS  I  use  R  AV  Cornell  's  nonlinear  compliance  model 
I  ref.  10)  that  formulates  tooth  stillness  as  a  function  of  position 
along  the  line  of  action.  This  compliance  model  is  based  on 
a  combination  of  the  stiffness  of  the  tooth  as  a  cantilever  beam, 
local  hertz  contact  compression,  and  fillet  and  tooth  foundation 
flexibility  effects  All  ot  the  above  except  the  local  contact 


compression  are  linear  functions  ot  the  load.  I  he  nonlinearity 
of  the  compliance  equation  is  due  to  the  hcrt/ian  deflection. 
PGT  uses  a  variable-variable  ,  >esh  stiffness"  (VVMS)  model 
tor  the  tooth  stiffness  The  VVMS  mode!  is  also  nonlinear  due 
to  local  hertz  contact  compression.  The  model  includes  tooth 
bending  effects  and  tooth  profile  errors  as  a  function  of  contact 
position. 

Tooth  root  stress  models.  — Of  the  four  programs 
investigated,  only  GRDYNMF'FT  and  DANST  are  capable 
ot  calculating  tooth  root  stresses.  Both  programs  use  the 
modified  Hevuood  formula  for  tooth  stress  sensitivity  as  given 
in  reference  1 1.  The  modified  Hey  wood  formula  calculates 
th  maximum  root  stress  as  a  function  of  tooth  contact  position, 
mesh  load,  face  w  idth,  stress  concentration  factor  of  the  fillet, 
and  basic  tooth  geometry.  The  formula  is  also  capable  of 
predicting  the  location  of  the  maximum  root  stress  on  the  tooth 
fillet.  The  modified  Heywood  formula  expresses  the  rixit  stress 
as  a  linear  function  of  the  applied  load  It  was  found  that  the 
formula  predicts  the  maximum  tensile  root  stress  w  ithin  about 
5  percent  of  finite-element  and  other  analysis  methods 
(ref  1 1 ). 

Input  error  models.— Actual  gear  systems  inherently  have 
one  or  more  types  of  errors  presen!.  In  an  attempt  to  more 
accurately  model  actual  systems,  all  of  the  programs  have 
prov  ided  some  means  of  including  errors  inherent  in  these 
systems.  The  program  GRDYNMUFT  allows  three  types  of 
errors  to  be  input.  These  are:  sun  runout  error,  helix  angle 
errors,  and  tooth  errors.  The  sun  runout  error,  applicable  to 
a  single  mesh  sy  stem  only  .  conv  erts  a  sun  center  displacement 
input  into  a  sinusoidal  tooth  spacing  error  array  to  simulate 
errors  associated  with  eccentrically  manufactured  gears.  The 
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ficurc  II  Program  DAVSI  riHslsI 

glc  errors  allow  the  user  to  input  a  constant  angular 
r  each  mesh  for  single  and  double  helical  gears.  The 
rors  are  comprised  of  tooth  error  arrays  on  five  teeth 
i  sun/planet,  ring/planet  mesh  This  tooth  error  input 
its  the  statistical  sunt  of  tooth  pitch  error,  profile  error, 
i  (or  planet  phasing)  error.  The  tooth  error  is  constant 
ic  profile  of  the  tooth.  The  program  POT  indirectly 


allows  two  types  of  error  to  be  input,  phase  error  and  tooth 
error.  The  phase  error  is  a  constant  lead,  or  lag.  tangential 
positioning  error  of  the  planets,  representing  planet  assembly 
inaccuracies.  The  tooth  error  consists  of  a  sinusoidal  error 
imposed  on  the  tooth  profile  with  the  amplitude  defined  by 
the  user.  This  error  models  gear  tooth  profile  manufacturing 
process  errors.  I  he  single  mesh  programs  (TELSGF.  and 
DANST)  have  available  tooth  profile  deviation  arrays. 
Deviations  from  the  true  involute  profile  can  be  defined  by 
inputting  the  corresponding  array  .  Tooth  spacing  error  can 
be  simulated  by  inputting  a  constant  deviation  along  the  tooth 
profile. 

Profile  modification.  —  Profile  modifications  are  often  used 
in  gears  to  lessen  engagement  impacts  in  attempts  to  reduce 
noise  and  vibration  in  gear  systems.  The  programs 
GRDYNMULT.  TELSGF.  and  DANST  allow  some  form  of 
modification  to  the  tooth  profile.  GRDYNMUI.T  incorporates 
an  equation  that  allows  the  user  to  input  the  deviation 
magnitude  at  the  tip.  length  of  the  modification  on  the  tooth 
profile,  and  the  shape  of  the  modification  curve.  To  determine 
the  profile  modification  curve  a  shape  factor  i,  input.  The 
default  shape  factor  (0)  produces  a  parabolic  profile 
modification.  A  linear  profile  modification  can  be 
approximated  with  this  equation  w  ith  a  shape  factor  of  -0.5. 
Other  shapes  associated  with  different  shape  factors  are  given 
in  reference  II  DANST  allows  two  standard  profile 
modifications  and  a  user  defined  shape  to  be  input.  A  standard 
linear  or  parabolic  tooth  profile  modification  can  be  chosen 
with  the  tip  dev  iation  magnitude  and  modification  length  along 
the  tooth  profile  input  by  the  user.  By  v  irtue  of  the  tooth  profile 
deviation  arrays  discussed  earlier,  other  user  defined  profile 


uodillcalions  can  ho  input  in  DAN^T.  P-  cram  TflSGT  also 
Hows  profile  modifications  to  ho  input  in  virtue  of  its  tooth 
notilo  deviation  array.  Standard  profile  modifications  such 
is  linear  and  parabolic  must  he  added  point  In  point  m  the 
irray 

Comparison  Runs  Study 

Short  of  using  experimental  data,  the  most  effective  wav 
>1  comparing  computer  programs  is  to  compare  their  output 
esults  based  on  common  input  values.  In  mis  study  the 
■>n 'grams  were  operated  using  common  model"  and  input 
xirameters.  Wheie  possible,  runs  where  performed  with 
parameters  such  as  speed,  load,  and  mesh  damping  vat ied  in 
ariu  r  to  obtain  program  comparisons  over  a  broad  spectrum 
.'I  input  coi'.  htions.  Input  parameters  eommon  to  at  least  two 
programs,  such  as  sun  center  stiffness,  were  also  varied  f<  r 
the  compart  on.  Due  to  the  nature  ot  the  programs,  two  types 
ol  input  models  were  required'  a  planet  try  system  with  three 
planets,  and  a  single  mesh  svstem.  A  discussion  of  the 
comparison  studv  results  are  thus  grouped  under  those  two 
categories 

Planetar)  Svstem  Runs 

Because  ot  the  svstem  limitations  of  the  program  PGT.  a 
three  plane!  planetar)  system  was  used  to  compare  programs 
P( i  I  and  GRDYNMl'I  .T.  I  able  VI  gives  a  description  of  the 
planetary  model  used  in  the  analysis.  along  with  the  undamped 
natural  frequencies  of  the  svstem.  as  calculated  In 
GRDYNMl'I.  I  As  seen  in  table  VI,  to  minimi/e  the  inlluenee 
of  the  input  and  output  peripherals  ol  PGT  in  the  analvsis. 
external  shaft  damping  and  mass  moments  of  inertia  were 
minimized,  and  external  shaft  stillness  values  were  maximized. 
Table  V  documents  the  comparison  runs  matrix  used, 
illustrating  which  parameters  were  varied  .u:j  their 
corresponding  values  Due  to  diliieulties  experienced  with  the 
p.ogram  PGT  and  with  the  HP  1000  computer  svstem.  only 
nine  comparison  runs  were  achieved.  Unfortunately  this  does 
not  allow  a  detailed  comparison  to  be  made:  however,  some 
general  observations  can  be  drawn.  Discussions  on  the  various 
parametric  rims  are  given  below. 

Speed  variation  runs. — To  compare  the  effect  of  input  speed 
on  the  maximum  dynamic  load  factor,  the  programs  were  run 
over  a  range  of  speeds  from  4000  to  8000  rpm.  figure  12  is 
a  plot  of  the  maximum  dynamic  load  factor  for  the  sun  planet 
mesh  .is  a  function  ol  input  speed,  as  predicted  by  both 
programs.  As  seen  in  this  figure,  both  programs  show  good 
correlation  except  at  GOOD  rpm  input  speed,  where  PGT 
predi.ts  a  peak  w  namic  load.  GRDYNMIT.T  preo  t'  a  peak 
at  the  2000  rpm  input  speed  point.  As  seen  in  table  VI.  this 
point  1 2000  rpm.  I 633  Hzi  is  within  7  percent  of  the  second 
harmonic  ot  the  second  natural  frequency  (1530  H/t.  as 
predicted  by  GRDYNMH.T  The  difference  in  pps'ieted  peak 


load  speeds  between  PG  T  and  GRDYN.VU'LT  could  be  due 
to  die  different  mesh  stillness  model  used  ;n  each  program, 
figure  17  is  the  same  plot  as  figure  12  except  that  the 
ma'i.mm’  ring  planet  mesh  loads  are  plotted.  Comparison  of 
figures  12  and  13  show  the  same  trends,  with  the  exception 
dial  the  ring  planet  plots  show  a  much  poorer  correlation 
between  the  two  programs. 

A  comparison  of  the  dynamic  mesh  load  p'ots  from  each 
program  through  one  tooth  mesh  cycle  at  input  speeds  of  4(XK). 
6000.  and  8()0t)  rpm  are  illustrated  in  figures  14.  15.  and  16. 
respectively  .  As  seen  in  figures  14  and  16.  the  sun  planet  mesh 
load  plots  are  very  similar  in  form  between  the  two  programs. 
The  ring  planet  mesh  load  plots  are  dissimilar  in  both  form 
and  magnitude,  f  igure  15  further  illustrates  the  discrepancy 
between  the  two  programs  at  the  6000-rpm  input  speed.  Here 
PGT  is  shown  to  predict  tooth  separation  with  a  maximum 
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dynamic  loud  factor  in  execs'-  of  2.8.  Again,  the  apparent 
difference  in  system  critical  speeds  could  he  due  to  different 
mesh  stiffness  models.  It  is  not  known  at  this  time  why  the 
ring  planet  mesh  loads  experienced  a  poorer  correlation  than 
the  sun  planet  mesh  loads. 

Sun  center  stiffness  runs. — The  sun  center  stiffness  input 
was  varied  in  each  program  to  compare  sun  center  flexibility 
effects  on  the  maximum  dynamic  load  factor.  Figure  17  plots 
the  relative  effects  on  the  maximum  dynamic  load  for  the 
sun  planet  mesh  for  three  sun  center  stiffness  values.  Three 
points  are  not  enough  to  provide  a  thorough  comparison; 


however,  some  general  trends  can  be  deduced  and  compared 
using  these  plots.  As  seen  in  figure  17.  trend  results  from  the 
two  programs  do  not  fully  agree.  PGT  favors  a  relatively  stiff 
sun  center  for  a  minimum  dynamic  load  factor,  whereas 
GRDYNMULT  indicates  an  optimum  sun  center  stiffness 
exists  between  the  two  extremes.  Similar  plots  for  the 
ring/planet  mesh  are  illustrated  in  figure  18.  Some  trends  can 
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Figure  17  — Comparison  ol  programs  PGT  and  ORDYNMIH.T  Maximum 
dwiarmc  load  (actor  as  a  function  of  sun  center  stiff  ness  for  the  sun/planet 
mesh  Input  torque  -  4.VV  N*ni:  input  speed  -  40<K)  rpm  (Table  V,  runs 
I .  h.  and  Vi 
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figure  }X  Comparison  ot  programs  PCiT  and  GKDYNMIT.T  Maximum 
ils  n.iiiiK  load  las  tor  as  a  function  o|  sun  center  stillness  lor  the  ring  planet 
mesh  Input  loiquc  '4  ‘l  Non.  input  sf  ed  4000  rpm  (Table  V.  runs 
I .  f>.  and  '0 


be  seen  in  this  figure:  however,  they  are  not  prominent  enough 
to  draw  any  conclusions. 

Damping  runs. —  The  mesh  damping  ratio  and  sun  center 
damping  coefficient  were  changed  to  compare  the  resulting 
effects  on  the  maximum  dynamic  load  factor  calculated  by  each 
program.  Figure  Id  illustrates  the  effects  on  the  maximum 
dynamic  load  factor  of  the  sun  planet  mesh  at  an  input  speed  of 
4tKK)  rpm  as  mesh  and  sun  center  damping  were  changed.  As 
seen  in  this  figure,  both  programs  show  an  increase  in  dynamic 
load  (9.0  percent  for  GRDYNMUL  f.  12.1  percent  for  PGT) 
as  the  mesh  damping  ratio  value  is  decreased  from  0.10  to 
0.03.  No  significant  change  was  noted  in  either  program  as 
the  sun  center  damping  coefficient  was  changed.  Similar  plots 
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Figure  70  C  omparison  of  programs  P(i  l  and  (iRDYNMI  T  T  Maximum 
dvuanik  load  factor  at  several  damping  conditions  for  the  ring  planet  mesh 
Input  torque  V*  0  Y*m.  input  speed  4000  rpm  (Table  V.  runs. 
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lor  the  ring:  planet  mesh  are  illustrated  in  figure  20.  The  trends 
seen  in  figure  20  are  similar  to  those  noted  in  figure  19: 
however,  they  are  not  prominent  enough  to  draw  any 
conclusions  from  them. 

Sun  center  movement.—  Sun  center  movement  is  calculated 
by  program  PGT  only,  thus  no  comparison  can  he  made  with 
GRDYNMULT.  PGT  predictions  of  the  sun  center 
displacement,  however,  proved  interesting  and  are  discussed 
below  .  Figure  21  illustrates  the  sun  center  movement  for  one 
revolution  at  a  variety  of  input  speeds.  The  maximum 
displacement  of  the  sun  center  is  seen  to  occur  at  6000-rpm 
input  speed,  the  same  as  with  the  maximum  dynamic  load 
factor.  As  the  speed  increases,  the  sun  center  displacement 
approaches  a  pattern  resembling  shaft  whirl.  As  expected,  the 
sun  center  movement  decreases  with  increasing  sun  center 
stiffness  (see  fig.  22).  A  decrease  in  mesh  damping  (from  0. 10 
to  0.03)  results  in  an  increase  in  sun  center  displacement  of 
more  than  two  times,  as  shown  in  figure  23.  Also  illustrated 
in  this  figure,  a  change  in  the  sun  center  damping  coefficient 
had  no  effect  on  the  sun  center  displacement  at  this  input  speed. 

Single  Mesh  Runs 

Because  of  the  system  limitations  ofDANST  and  TELSGE. 
a  single  mesh  system  was  used  to  compare  programs 
GRDYNMULT.  TELSGE.  and  DANST.  Table  VII  gives  a 
description  of  the  single  mesh  model  used  in  the  analy  sis  along 
with  the  undamped  natural  frequencies  of  the  system  calculated 
by  each  program.  As  seen  in  this  table,  the  programs  predicted 
similar  natural  frequencies  for  the  single  mesh  system  (all 
w  ithin  13  percent  of  the  calculated  average  of  4532  Hz).  Of 
the  three  programs,  only  DANST  includes  external  shafts  and 
masses  in  the  system  dynamics.  To  maintain  an  equal 
comparison  basis  among  the  three  programs,  it  was  necessary 
to  minimize  the  influence  of  the  peripheral  masses  in  program 
DANST.  This  was  accomplished  by  using  highly  flexible  input 
and  output  shafts  in  the  program.  In  the  planetary  sy  stem  runs 
program  PGT  used  short,  highly  rigid  shafts  with  small 
peripheral  mass  inertias  to  minimize  their  effects  on  the  system 
dynamics.  This  method  did  not  work  as  well  with  program 
DANST.  thus  the  opposite  approach  of  flexible  shafts  was  used 
to  isolate  the  peripheral  mass  inertias  from  the  mesh  dynamics, 
figure  24  illustrates  the  effect  of  varying  the  magnitude  of 
the  peripheral  masses  on  the  maximum  dynamic  load  factor, 
as  predicted  by  program  DANST  with  the  flexible  shaft 
configuration.  As  seen  in  this  figure,  the  dy  namic  load  factor 
changes  minimally  with  peripheral  mass  changes,  indicating 
good  isolation  of  the  mesh  dy  namics  with  this  configuration. 
Table  VII)  documents  the  comparison  runs  matrix  used, 
illustrating  which  parameters  were  varied  and  their 
corresponding  values.  Discussions  comparing  the  effects  of 
the  various  parametric  runs  on  the  variables  calculated  by  the 
programs  are  given  below. 

Dynamic  load  factor. — A  v  ariety  of  input  speeds  and  torques 
were  used  to  compare  the  relative  effects  of  speed  and  load 
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igure  22  —  PGT  program  sun  tenter  movement  predictions  at  various  sun 
center  stiffness  values  Input  torque  =  33.9  N*m;  input  speed  =  4000  rpnt 
(Table  V.  runs  1.  6.  and  9). 
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figure  23  —  PGT  program  sun  center  movement  predictions  for  various 
damping  values  Input  torque  =  33.9  N«m;  input  speed  =  4(KX)  rpm 
(Table  V,  runs  I.  7.  and  X) 
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Figure  24  —Effect  on  maximum  dynamic  load  factor,  as  function  of  input 
speed,  as  peripheral  mass  multiplication  factor  (PMMF)  is  varied  in  program 
DANST,  with  highly  flexible  input  and  output  shafts.  < Jlnrul  ma„.  =  PMMF 
*^Aming  star'  Lulpul  mj"  —  PMMFXjjmcn  pt..,r  I . 


Figure  25. -Comparison  of  program  GRDYNMULT.  DANST.  and  TELSGE. 
Maximum  dynamic  load  factor  as  function  of  input  speed.  Input  torque 
=  203.4  N*m  (Table  VIII,  run  I). 


on  the  dynamic  load  factor  as  calculated  by  each  program. 
Maximum  dynamic  load  factors  are  plotted  as  a  function  of 
input  speed  for  an  input  load  of  203.4  N»m  (1800  in.«!b)  in 
figure  25.  As  seen  in  this  figure,  all  three  programs  show  good 
correlation  (average  difference  within  5  percent)  except  at 
5500  rpm,  where  TELSGE  results  diverge.  This  speed  is 
within  8  percent  of  the  speed  corresponding  to  the  half 
harmonic  of  the  natural  frequency  predicted  by  TELSGE 
(5130  rpm).  This  half  harmonic  phenomenon  is  also  seen  in 
programs  GRDYNMULT  and  DANST,  although  at  a  lesser 
degree.  DANST  and  GRDYNMULT  both  indicate  peaks  at 
the  5000-rpm  data  point.  The  predicted  half  harmonic  speed 
of  program  DANST  (5191  rpm)  is  within  4  percent  of  this 
peak  dynamic  load  point.  The  corresponding  half  harmonic 
speed  of  program  GRDYNMULT  (4246  rpm)  is  within  15 
percent  of  the  peak  dynamic  load  point.  Because  the  mesh 
stiffness  varies  with  tooth  position  during  mesh,  the  predicted 
natural  frequencies  are  only  estimates  of  the  actual  values, 
based  on  assumed  constant  mesh  stiffness  quantities.  A 
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Figure  26.— Comparison  of  programs  GRDYNMULT.  DANST.  and  TELSGE. 
Dynamic  load  factor  as  function  of  contact  position.  Input  torque 
=  203.4  N*m:  input  speed  =  2000  rpm  (Table  VIII.  run  5). 
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Figure  27.—  Comparison  of  programs  GRDYNMULT.  DANST,  and  TELSGF. 
Dynamic  load  factor  as  function  of  contact  position.  Input  torque 
=  203.4  N»m;  input  speed  =  4000  rpm  (Table  VIII,  run  7). 
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Figure  28— Comparison  of  programs  GRDYNMULT.  DANST.  and  TF.LSGE. 
Dynamic  load  factor  as  function  of  contact  position.  Input  torque 
=  203.4  N«m;  input  speed  =  6000  rpm  (Table  VIII.  run  9). 
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Figure  29  —Comparison  of  programs  GRDYNMULT.  DANST,  and  TELSGE. 
Dynamic  load  factor  as  function  of  contact  position.  Input  torque 
=  203  4  N*m;  input  speed  =  8000  rpm  (Table  VIII,  run  1 1) 
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Figure  30.— Comparison  of  programs  GRDYNMULT.  DANST.  and  TELSGE. 
Maximum  dynamic  load  factor  as  function  of  input  torque  at  2000-rpm 
input  speed  (Table  VIII.  runs  5.  17-20). 
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Figure  31— Comparison  of  programs  GRDYNMULT.  DANST.  and  TELSGE. 
Maximum  dynamic  load  factor  as  function  of  input  torque  at  6000-rpm 
input  speed  (Table  VIII.  runs  12-16). 

comparison  of  the  actual  dynamic  load  plots  from  each 
program  for  a  variety  of  speeds  can  be  seen  in  figures  26  to 
29.  As  illustrated  in  these  figures,  the  dynamic  load  factor  plots 
are  very  similar  in  both  magnitude  and  form.  All  three 
programs  show  a  decrease  in  the  frequency  of  dynamic  load 
fluctuations  as  the  input  speed  increases,  and  a  condition  close 
to  tooth  separation  at  the  8000  rpm  input  speed  (fig.  29). 
Figures  30  and  3 1  are  plots  of  the  maximum  dynamic  load 
factor  as  a  function  of  input  torque  for  input  speeds  of  2000 
and  6000  rpm.  respectively.  As  seen  in  these  figures,  the 
programs  predict  a  fairly  constant  dynamic  load  factor  regard¬ 
less  of  the  input  torque  value.  This  is  as  expected  since  the 
dynamic  and  static  load  are  both  linear  functions  of  the  input 
torque. 

Tooth  root  stress.— Tooth  root  stress  was  another  variable 
compared  using  a  variety  of  input  loads  and  torques.  As  illus¬ 
trated  in  figure  32,  the  maximum  root  stress  predicted  by 
DANST  and  GRDYNMULT  correlate  reasonably  well  through 
the  speed  range,  showing  similar  form  and  magnitudes  that 
disagree  only  slightly  (average  difference  within  16  percent). 
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Figure  32.—  Comparison  of  programs  GRDYNMUI.T  and  DANST.  Maxi¬ 
mum  loofh  root  stress  as  function  of  input  speed.  Input  torque 
=  203.4  N*m  (Table  VIII,  runs  5-11). 
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Figure  33  —Comparison  of  programs  GRDYNMULT  and  DANST.  Maxi¬ 
mum  tooth  root  stress  as  function  of  input  torque  at  2(XX)-rpm  input 
speed  (Table  VIII.  runs  5,  17-20). 


As  expected,  both  show  peak  values  at  the  5000  rpm  data  point. 
Figures  33  and  34  plot  the  maximum  tooth  root  stress  as  a 
function  of  input  torque  at  input  speeds  of  2000  and  6000  rpm. 
respectively.  As  seen  in  these  figures,  both  programs  show 
the  tooth  root  stress  to  be  relatively  linear  with  input  torque. 
This  is  expected  since  both  use  a  form  of  the  modified 
Heywood  formula  which  gives  tooth  root  stress  as  a  linear 
function  of  applied  load. 

Contact  stress.— The  local  contact  pressure,  or  hertz  stress, 
is  calculated  by  programs  TELSGE  and  GRDYNMULT.  As 
seen  in  figure  35,  both  programs  show  similar  trends  and 
values  (average  difference  within  4  percent)  with  input  speed 
with  the  exception  of  the  TELSGE  results  between  5000  and 
6000  rpm.  Here,  due  to  the  close  proximity  of  the  half 
harmonic  of  the  system.  TELSGE  would  not  converge.  Both 
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Figure  34  — Comparison  of  programs  GRDYNMUI.T  and  DANST.  Maxi¬ 
mum  tooth  root  stress  as  function  of  input  torque  at  6000-rpm  input 
speed  (Table  VIII.  runs  0.  12-lb). 
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Figure  35.— Comparison  of  programs  GRDYNMULT  and  TELSGE.  Maxi¬ 
mum  Hertz  stress  as  function  of  input  speed.  Input  torque  =  203.4  N«m 
(Table  VIII.  runs  5- 1 1 ), 
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Figure  3b  — Comparison  of  programs  GRDYNMUI.T  and  TELSGE.  Maxi¬ 
mum  Hertz  stress  function  of  input  torque  at  2000-rpnr  input  speed 
(Table  VIII.  runs  5.  17-20). 
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Figure  37.  -Comparison  of  programs  GRDYNMUI.T  and  THLSGE.  Maxi¬ 
mum  Hertz  stress  as  funetion  of  input  torque  al  6000-rpm  input  speed 
(Table  VIII.  runs  12- 16). 
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Figure  38.— Comparison  of  programs  GRDYNMULT  and  TELSGE.  Maxi¬ 
mum  Hash  temperature  as  function  of  input  speed.  Input  torque 
=  203.4  N.m  (Table  VIII.  runs,  5-11). 


programs  predicted  nearly  identical  trends  and  values  with 
input  torque  variations,  as  seen  in  figure  36  for  a  2000-rpm 
input  speed  and  figure  37  for  a  6000-rpm  input  speed.  The 
nonlinear  relationship  between  input  torque  and  hertz  stress 
can  be  clearly  seen  in  figures  36  and  37. 

Flash  temperature.— The  flash  temperature,  as  calculated 
by  programs  TELSGE  and  GRDYNMULT,  was  the  last 
variable  compared  using  a  variety  of  input  torques  and  speeds. 
Generally,  it  was  found  that  both  programs  predicted  similar 
trends  with  input  speed  and  input  torque;  however,  actual 
values  differed  by  between  46  and  153  K  (83  and  275  °F). 
Figure  38  illustrates  the  similar  speed  trends  displayed  by  both 
programs.  TELSGE  did  not  converge  in  the  input  speed  region 
between  5000  and  6000  rpm.  Maximum  flash  temperatures 
are  plotted  as  a  function  of  input  torque  in  figures  39  and  40 
at  2000-  and  6000-rpm  input  speeds,  respectively.  As  seen 
in  these  figures,  both  programs  displayed  the  same  nonlinear 
increasing  flash  temperature  trend  with  increasing  input  torque. 

Profile  modification.— To  compare  the  relative  effects  of 
profile  modification  on  the  dynamic  load  factor  as  calculated 
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Figure  39.  -Comparison  of  programs  GRDYNMULT  and  TELSGE.  Maxi¬ 
mum  Hash  temperature  as  function  of  input  torque  at  2(KK)-rpm  input 
speed  (Tab)e  VIII.  runs  5.  17-20). 
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Figure  40. —Comparison  of  programs  GRDYNMULT  and  TELSGE.  Maxi¬ 
mum  Hash  temperature  as  function  of  input  torque  at  6000-rpm  input 
speed  (Table  VIII.  runs  12-16). 


Figure  4 1 .  —Comparison  of  programs  GRDYNMULT.  DANST.  and  TELSGE 
Maximum  dynamic  load  factor  as  function  input  speed,  with  tooth  profile 
modification  Input  torque  =  203.4  N»m  (Table  VIII.  run  3). 

by  each  program,  a  standard  tip  relief  was  added  to  the  single 
mesh  system.  The  tip  relief  consisted  of  a  parabolic  shape  along 
50  percent  of  the  length  from  the  tip  to  the  pitch  point,  with 
a  maximum  deviation  magnitude  of  0.0178  mm  (0.0007  in.) 
at  the  tip.  Plots  of  the  dynamic  load  factor,  as  a  function  of 
input  speed,  with  profile  modification  are  given  in  figure  41 . 
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Comparison  of  figure  4!  w  ith  figure  25  (same  run  parameters 
as  fig.  41  hut  with  no  profile  modification)  shows  that  the  most 
dramatic  amplitude  reductions  occur  similarly  in  programs 
DANST  and  TELSGE  at  speeds  near  their  predicted  half  har¬ 
monic  speeds  DANST  shows  an  amplitude  reduction  of  35 
percent  at  the  5000-rpm  data  point  (predicted  half  harmonic 
speed  at  5191  rpm).  TELSGE  reduces  from  a  divergance 
situation  to  a  maximum  dynamic  load  factor  of  1 .27  with  profile 
modification  at  the  55(X)-rpm  data  point  (predicted  half  harmonic 
speed  at  5130  rpm).  TELSGE  and  DANST  also  experienced 
similar  dynamic  load  factor  reductions  at  speeds  below  the  peak 
amplitude  speeds  with  profile  modification  added,  as  illustrated 
in  tigure  41  GRDYNMULT  showed  no  appreciable  difference 
with  profile  modification  added.  It  is  not  known  at  this  time 
why  GRDYNMIT.T  did  not  show  any  change  with  the  addition 
of  profile  modification  in  this  example. 

Mesh  damping.—  To  compare  the  relative  effects  of  the 
mesh  damping  ratio  on  the  dynamic  load  factor,  a  number  of 
runs  were  made  with  mesh  damping  ratio  input  values  ranging 
from  0.03  to  0. 17.  Because  damping  effects  are  more  promi¬ 
nent  at  system  resonanee'points,  an  input  speed  of  5000  rpm 
was  chosen  because  of  its  close  proximity  to  the  half  harmonic 
speeds  predicted  by  each  program.  As  illustrated  in  figure  42. 
all  three  programs  show  good  correlation  at  damping  ratios 
of  0. 10  or  greater.  As  seen  in  this  figure,  all  of  the  programs 
predict  a  reduction  in  maximum  dynamic  load  factor  as  the 
mesh  damping  ratio  value  is  increased  from  0. 10  to  0.17  <  12 
percent  reduction  for  TELSGE.  19  percent  reduction  for 
GRDYNMU LT.  and  14  percent  reduction  for  DANST).  At 
damping  ratios  lower  than  10  percent,  the  TELSGE  program 
diverged.  The  close  proximity  of  the  5000-rpm  input  speed 
to  the  hall  harmonic  of  the  first  natural  frequency  predicted 
by  TELSGE  (within  3  percent  of  5 1 30  rpm)  is  most  probably 
the  reason  TELSGE  is  highly  sensitive  to  the  mesh  damping 
ratio  changes  at  this  speed.  DANST  and  GRDYNMULT  show 
good  correlation  over  the  w  hole  range  of  damping  ratios  used. 
As  seen  in  figure  42.  as  the  mesh  damping  ratio  increases  from 
0.03  to  0.17  both  programs  show  a  near  identical  decrease 
of  the  dynamic  load  factor  in  both  form  and  magnitude 
(DANST:  22  percent  reduction.  GRDYNMULT:  23  percent 
reduction)  from  an  average  value  ol  1.64  to  1.27. 
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Concluding  Remarks 

A  comparison  study  was  performed  w  ith  the  gear  dynamic 
analysis  computer  programs  PGT.  GRDYNMULT.  TELSGE, 
and  DANST  at  NASA  Lewis  Research  Center  The 
comparison  study  consisted  of  two  major  parts.  The  first  part 
inv  olved  a  direct  comparison  of  the  capabilities,  input  options, 
and  output  options  of  the  programs.  Results  of  this  study  w  ere 
tabulated  and  some  general  comments  are  as  follows: 

1.  GRDYNMULT  appears  to  be  the  most  versatile  in 
system  si/e.  ty  pe,  and  analysis  capabilities  of  all  the  programs 
compared. 

2.  TELSGE  provides  the  most  detailed  analysis  on 
lubrication  dynamics,  y  ielding  quantities  such  as  film  thickness 
and  Hash  temperatures. 

3.  DANST  incorporates  the  most  versatile  tooth  profile 
deviation  routine,  allow  ing  the  user  to  enter  standard  or  user 
defined  shapes  and  magnitudes. 

4.  PGT  provides  a  sun  center  movement  routine  which 
allow  s  the  user  to  obtain  the  displacement  of  the  sun  center 
through  one  or  more  revolutions. 

The  second  part  of  the  comparison  study  involved 
performing  parametric  comparison  runs  using  identical  input 
models.  Some  general  results  from  this  study  are  given  below : 

1 .  Computer  programs  PGT  and  GRDYNMULT  predicted 
similar  levels  and  form  of  the  dynamic  sun/planet  mesh  loads 
as  the  input  speed  was  varied.  Ring/planet  mesh  loads  differed 
significantly  between  the  programs. 

2.  Programs  TELSGE.  GRDYNMULT.  and  DANST  all 
predicted  dynamic  mesh  loads  of  similar  form  and  magnitudes 
as  the  input  speed  and  torque  were  varied.  TELSGE  results 
diverged  at  input  speeds  near  its  half  harmonic  resonant  speed. 

3.  Root  stress  predictions  from  programs  DANST  and 
GRDYNMULT  showed  good  trend  correlation  with  input 
speed  and  torque  variations.  Magnitudes  correlated  reasonably 
well  with  only  minor  variations. 

4.  Programs  TELSGE  and  GRDYNMULT  predicted  nearly 
identical  hertz  stress  levels  and  trends  us  input  torques  and 
speeds  were  varied. 

5.  Programs  TELSGE  and  GRDYNMULT  predicted 
similar  Hash  temperature  trends:  however,  actual  values  were 
not  in  close  agreement.  GRDYNMULT  consistently  predicted 
higher  than  expected  Hash  temperatures. 

Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland.  Ohio.  Decemher  14.  14X8 
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TABLE  I  -GENERAL  INFORMATION 


Program  name 

Dynamics  of 
planetary 
gear  trains. 
PGT 

Epicyclic  gear 
dynamic  analysis 
program. 
GRDYNMULT 

Thermal  elasto- 
hydrodynamtc 
lubrication  of 
spur  gears. 
TELSGE 

Dynamic  analysis 
of  spur  gear 
transmissions. 

DANST 

Documentation 

1 

NASA  CR  3793 

NASA  CR- 179563 

NASA  CR-3241. 
NASA  TP-2610 

NASA  CR- 179473 

Operating  system 

HP  1000 

IBM  370 

Cray  XMP 

IBM  370 

System  types 

Planetary  with 
input' output 
peripheral 
components 

Planetary  system,  star 
system,  differential 
system,  single  mesh 
(external-external), 
single  mesh 
(external-internal) 

Single  mesh 
(external-external) 

Single  mesh 
(external-external) 
with  input  and 
output  peripheral 
components 

Ssstem  si/e 

Single  stage 
planetary  system 
with  three  planets 

Single  mesh, 
epicyclic  gear  train 
with  20  planets 
maximum 

Single  mesh 

Single  mesh 

Gear  types 

Spur  gears 

Spur  gears,  single 
helical,  double 
helical 

Spur  gear 

Spur  gear 

Tooth  forms 

- - - 

Internal,  external 

Internal,  external, 
buttress 

External 

External 

Maximum  contact 

ratio 

2.0 

3.0 

2.0 

2.0 

TABLE  II. -COMPARISON 


Program 

Gear  life 
calculation 

Calculates 
dynamic  mesh 
load 

Calculates 

combined 

stiffness 

!  Calculates 

tooth  root 

stress 

! 

Calculaies 

moth  hertz 

stress 

Film 

thickness 

calculations 

Calculates 

surface 

temperature 

PGT 

No 

Yes.  tor  each  mesh 
at  each  planet,  as  a 
function  of  roil  angle 

i - 

Yes.  tor  each  mesh 
at  each  planet,  as  a 

I  function  of  roll  angle 

No 

No 

No 

No 

ORDYNMLLT 

No 

Yes.  maximum  value 
for  each  mesh  at 
each  planet,  and  as 
a  function  of  position 
along  line  of  contact 

!  Yes.  compliance 
i  function  coefficients 

1  calculated 

l 

-i - 

Yes.  maximum 
j  value  tor  each 
|  gear  in  mesh, 
and  as  a  function 
of  position  along 
line  of  contact 

Yes.  maximum 

value  for  each 

mesh  at  each 
'  planet,  and  as 
a  function  ot 
position  along 
line  of  contact 

No 

No 

TELSGE 

Yes.  based 
on  dynamic 
;  mesh  loads 

Yes,  maximum  value. 

1  and  ax  a  function  of 
i  position  along  line  ot 
!  contact 

|  Yes.  as  a  function  ol 

1  position  along  line  ot 
contact 

'  No 

! 

1 

l 

i 

Yes.  as  a 
function  of 
position  along 
line  of  contact 

Yes.  as  a 

lunction  of 
position  along 
line  of  contact 

Yes.  gear  and 
pinion  surface 
temperature  *s 
a  function  ot 
position  along 
line  of  contact 

D  AN  ST 

No 

Yes.  as  a  function  of 
roll  angle 

Yes.  as  a  function  of 
roll  angle 

Yes .  tor  each 

1  gear  as  a  function 
roll  angle 

No 

No 

No 

>F  CAi’ABII.ITU.S 


Calculations 

flash 

»cmpe  rat  u  re 

Parameter 

run  survev 

Calculates 
dynamic  PV 
surface  damage) 

factor 

Geometric 

preprocessor 

Natural 

frequency 

predictions 

frequency 

analysis 

Sun  center  movement' 

calculation 

No 

No 

No 

Yes.  tor  tooth  geometry 
input  in  dynamic  load 
calculation 

No 

No 

Yes 

Yes.  maximum 

v alut  tor  ejeh 

mesh  at  each 
planet,  and  as  a 

funct  ton  ot 

poMt  ion  alone 
line  •  >t  contact 

Yes,  speed  run. 
determines  maximum 
dynamic  load  jt  each 
mesh  tor  eacn  speed 
increment 

Yes.  maximum 

v  due  for  each 

mesh  at  each 
planet,  and  as  a 
function  ot 
position  along 
line  of  contact 

Yes.  lor  ,.>oth  geometry 
input  in  dynamic  load 
calculation,  and 
determines  optimum 
profile  modification 

Yes.  predicts 
natural  frequencies  1 
jnd  mode  shapes  ot  i 
the  system 

Yes.  with  post¬ 
processing  program 
"I replot*'  performs  i 
trequency  analysis 
on  mesh  load 

v  anations 

No 

Yes.  gear  id 
pint*)  n  Hj^h 
•em^erature  as  a 

‘unci  Jon  ot 

•s  «mi  ion  atone 

ne  •  »f  contact 

Yes.  speed  run 
load  run.  face  Aidth 

run.  outside  radius 

run.  number  ot  teeth 

run.  and  surface 

convection  heat 

transfer  coefficient  run 

No 

Yes.  for  looth  geometry 
input  in  dynamic  load 
calculation 

Y  es.  predicts  system  1 
natural  trequency 

No  j 

1 

N/A 

No 

Y-.'v  -peed  run. 
determines  maximum 

dsnumis  load  .it  each 

-peed  increment 

No 

Yes.  tor  tc  th  geometry 
input  in  dynamic  load 
calculation 

Yes.  predicts 
first  three  natural  i 
frequencies  ot  the 
system 

Yes.  performs 
frequency  analysis 
on  dynamic  mesh 
load,  and  on  the 

static  transmission 

error 

N  A 

TABLE  III. -PROGRAM 


Program 

lubrication 

Iteration 

;  convergence 
;  tolerance 

Mesh  damping 
ratio 

Planet  gears 
phasing 

constant 

Face  width 
crowning 
parameter 

Stiffness 
of  peripheral 
shafts 

Damping  1 

of 

peripheral  ) 
shafts  | 

PGT 

No 

No 

Yes 

| 

; 

No 

No 

Yes.  actual 
stiffness 
value  entered 

Yes.  damping 
coefficient 

entered 

GRDYNMll.T 

Ves.  choice  of 
several  oils  in 
program,  or 
user  defined  oil 

Yes.  number  of 
i  iterations,  and 

convergence 
,  tolerance 
!  can  be  input 

Yes 

Yes 

Yes.  length  of 
face  width  crown, 
and  edge  relief 
are  input 

No 

No 

ILLSGh 

Yes.  user  inputs 
oil  tspe  and 
properties 

1 - 

1  No 

Yes 

1 

N  A 

No 

No 

No 

DANST 

Ves.  user  can 
define  one  ol  two 
lubrication  models 

available 

No 

Yes 

i 

! 

N'A 

No 

Yes,  user  inputs 
shaft  diameter, 
length,  and 
modulus  for  both 
;  input  and  output 
shafts 

No 

IN  PC  r  OPTIONS 


|  .oa«J  and  Gear 

di  tver  mass  material 

foment  of 

inertia 

- . - 

Errors  Profile 

modifications 

Planet 

carrier 

flexibility 

Ring  gear  Floating  sun 

rim  flexibility  gear 

: 

Vos  Yes.  Young's  modulus. 

Poisson's  ratio,  material 
density 

Y  es.  planet  phase  1  None  indicated  in 
angle  error  input.  1  documentation 
tooth  profile  error  j 
input 

_ 1 _ 

No 

Yes.  sun  center 
stiffness  and 
damping  coefficient 
can  be  input 

No  Yes.  Y  oung's  modulus. 

Poisson's  ratio,  material 
density 

Yes,  tooth  spacing 
errors  input,  sun 
gear  run-out 
errror.  heli  angle 
error 

Yes.  tooth  profile 
modification  shape, 
length  along 
tooth  surface,  and 
magitude  can  be 
input 

Yes.  azimuthal 
planet  carrier 
stiffness  can  be 
input 

Yes,  azimuthal 
ring  gear  nm 
stiffness  can  be 
input 

Yes.  sun  center 

stiffness  and 
damping  coefficient 
can  be  input 

No  Y  es.  Young's  modulus. 

|  Poisson  s  ratio,  material 
|  density,  specific  heat. 

!  thermal  conductivity 

Yes.  100  point 
array  available  for 
user  defined  tooth 
profile  deviation 

Y'es.  100  point 
array  available  for 
user  defined  tooth 
profile  deviation 

N/A 

N'A  N/A 

| 

; 

i 

Y  es  Yes.  Young's  modulus 

1  Poisson's  ratio,  material 
|  density 

Yes,  121  point 
array  available  for 
user  defined  tooth 
profile  deviation 

Yes.  tip  relief 
parameters  can  be 
input  for  a  linear 
or  parabolic 
shape,  or.  user 
can  define  shape 
using  121  point 
file  deviation  array 

N/A 

N/A  N/A 

:j 


1  ABLF.  IV  -PROGRAM  Ol  TPl  T  OPTIONS 


Pros:  ram 

POT 

GRDVNMl  l.r 

IFI.SGF 

DA. VST 

Printed 

Printed  as  a  luncti. 

*n  ot 

Geometric  preprocessor  results 

Input  data  Hcho 

Input  data  hcho 

output 

gear  mil  angle 

Involute  modification  tables 

Ciear  lite  calculations 

Gear  teeth  deflection 

i\  j  liable 

•  Mesh  stillness 

Input  data  Hcho 

Printed  as  a  tunction  ot 

Static  transmission  error 

•  Vlc>h  dMUIllk  . 

•  \lds 

Constants  tor  the  lourth  order 

c<  intact  position 

Dxnainic  tooth  load 

compliance  1  unction 

•  Combined  stittness 

Fourier  transform  ot  the  dvnamtc 

Boundarv  conditions  iteration 

•  Dvnamtc  load  t actor 

mesh  loads 

results 

•  Hertz  stress 

Fourier  transform  ot  the  static 

Maximum  values  tor  each  mesh: 

•  Filin  thickness 

transmission  error 

•  Hertz  stress 

•  Pinion  temperature 

•  Root  stress 

•  Gear  temperature 

•  Dvnamtc  load  factor 

•  Flash  temperature— pinion 

•  Plash  temperature 

•  Flash  temperature— eear 

•  Dynamic  P\ 

Other  values  printed 

Sv  vie m  natural  trequenev  results 

•  Mesh  natural  trequenev 

•  Maximum  dvnamtc  load 

•  Average  mesh  stillness 

Plotted 

Plotted  as  a  lunette 

*n  M 

Plotted  as  a  1  unction  of  position 

Plotted  routine  not  available 

Plotted  as  a  tunction  of  gear  roll  angle. 

output 

jo  sir  roil  angle 

along  line  of  action 

on  Cray  version 

•  Tooth  dellectton 

available 

•  Mesh  stillness 

•  Dynamic  load  factor 

•  Static  transmission  error 

•  Mesh  dsnamic  . 

a  Mils 

•  Pressure  slidine  velocity  tPVi 

•  Tooth  stittness 

Sun  gear  center  tin 

h  etnent 

•  Hertz  stress 

•  Tooth  load  sharing 

•  Flash  temperature 

•  Coefficient  ot  friction 

•  Root  stress— each  eear 

•  Torsional  torque 

Frequency  analysis  ol  mesh  loads 

•  .Static  and  dynamic  tooth  lojds 

at  each  mesh 

•  Tooth  root  stress— gear  and  pinion 

l-'t'uricr  Irunsiorm  plot  ot  dynamic  mesh 
loads 

Fourier  transform  plol  ot  the  static 
transmission  error 
Dynamic  load  factor  plot  for  speed 
survev  run 


FAB1.H  V  — ORDYNM L'LT-PGT  COMPARISON  Rl'NS  MATRIX  INSCRIPTION 
[Input  torque  =  0  N»m  i.'OO  1  n •  i h >  lor  ail  runs  | 


Run 

number 

Input  speed, 
rpm 

Sun  center 

stiffness 

Mesh  damping 
ratio. 

Sun  center 
damping 

N  m 

lb  in. 

\«s  m  lb-s  m 

i 

4t  KKI 

l'51  0  ■  10" 

Ill  ■  10" 

o  1 

r  ,m  "i 

: 

5000 

bOOO 

4 

“000 

5 

soon 

t 

♦ 

b 

4000 

l~5  1  ■  10"  1 

I  O  ■  JO*' 

»  » 

*» 

P5I Ox  |o"  1 

10.  UP 

* 

y  50 

X 

l“5l  Ox  10" 

lo  «  or 

(M 

r  5i  i 

d 

*  1 

I'M  -  III" 

o  i  ■  !u‘ 

i 

1 '  5 1  I 

TABLE  VI. -SYSTEM  DESCRIPTION  AND  NATURAL 
FREQUENCY  PREDICTION  OF  THE  PLANETARY 
GEAR  TRAIN  USED  IN  PGT— GRDYNMULT 
COMPARISON  RUNS 


System  description 

Ssstem  tv  pc . 

Planetary 

Diametral  pitch . 

....  8.4667 

Pressure  angle,  deg ... 

22  5 

Number  of  teeth 

Sun  . 

.  14 

Planets  . 

. 28 

Ring  .  .  . 

. 70 

Number  of  planets  .... 

.  .  .  3 

Face  width. 

mm  (in.) 

Sun 

30  (1.181 1) 

planets  . 

30  (1.1811) 

ring 

36  (1.4184) 

Natural  frequency  predictions  (from  program  GRDYNMULT) 

N 

fn.  Hz 

2.fn,  Hz 

3«fn.  Hz 

4*fn.  Hz 

1 

144 

288 

432 

576 

*> 

765 

1530 

2295 

3060 

3 

1020 

2040 

3060 

— 

4 

1416 

2832 

— 

— 

5 

2378 

— 

— 

6 

25 1 3 

— 

— 

For  PGT  input  only 

J  (driven.  N*m*s;  (in.«lb«s-) . 113x10  6  (0.001) 

J  (load).  N*m*s-  (in.«lb»s-) . 1 13x  10^h  (0.001 ) 

Inpul  shaft  stiffness.  N*m  (lb-in.) .  l^.SxlO11  (10x10'’) 

I  Output  shaft  stiffness.  N»m  ( I b  Tn  ) .  I  75  x  1 0V  1 1 0  x  1 ) 

i  Input  shaft  damping.  N«sm  (lh»s,in.) .  0.175  (0.001 1 

j  Output  shaft  damping.  N«s  m  (lb»s'in)  .  0.175  (0.001 1 


TABLE  VII.— SYSTEM  DESCRIPTION  AND  NATURAL 
FREQUENCY  PREDICTIONS  OF  THE  SINGLE  MESH 
SYSTEM  USED  IN  GRDYNMULT— TELSGF.—DANST 
COMPARISON  RUNS 


System  description 


System  type . Single  mesh 

Diametral  pitch . 8.000 

Pressure  angle,  deg . 20 

Number  of  teeth  (pinion) . 28 

Number  of  teeth  (gear) .  28 

Face  width,  mm  (in.) . 6.35  (0.25) 

Lubrication . MIL-L-23699 


Natural  frequency  predictions 


1 

GRD'tNMULT 

TELSGE 

1 

w.  rpm 

fn.  Hz 

w,  rpm 

1 

fn.  Hz 

w.  rpm 

1 

71 

3963 

8492 

4788 

10  200 

2 

40 

86 

— 

— 

-  | 

— 

3 

4845 

10  382 

— 

— 

-  , 

— 

For  DANST  input  only: 

J  (driver) . I*J  (pinion) 

J  (l-std)  .  M  (gear) 

Input  shaft  diameter,  mm  (in.) .  5.08  (0.20) 

Output  shaft  diameter,  mm  (in  )  .  5.08  (0.20) 

Input  shaft  length,  mm  (in.)  . 381  (15.0) 

Output  shaft  length,  mm  (in.) . 381  (15,0) 
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TABLE  VIII— DANST-GRDYNMULT-TELSGE  COMPARISON 
RUNS  MATRIX  DESCRIPTION 
| A  =  GRDYNMULT.  B  =  DANST.  and  C  =  TELSGE  ] 
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